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Optical  excitation  of  electrons  in  an  asymmetric  double  quzmtum  well  is  theoretically 
examined.  The  well  is  biased  to  align  the  excited  levels  and  permit  resonant  electron  tunnel¬ 
ing.  Emphasis  is  made  on  the  photoinduced  transfer  of  electrons  counter  to  the  bias  electric 
field  force.  Systems  with  a  weak  polarization  relaxation  (dephasing)  are  studied  using  the  con¬ 
ventional  technique  of  the  Schrodinger  equation.  A  density  matrix  approach  is  developed  to 
describe  optical  excitations  in  the  presence  of  an  arbitrary  dephasing.  Quantum  beats,  which 
follow  a  short-pulse  excitation  of  the  double  well,  are  shown  to  crucially  depend  on  the  de¬ 
phasing.  The  excitation  profiles  obtained  for  cases  of  different  dephasing  reveeJ  the  full  range 
of  tunneling  coupling  between  the  wells  from  completely  coherent  to  incoherent  (stepwise). 


PACS:  73.40.Gk,  73.20.Dx,  73.50.Pz,  73.40.Kp 


INTRODUCTION 


The  aim  of  this  paper  is  to  theoretically  consider  processes  of  opticeil  excitation  and 
electron  transfer  in  an  asymmetric  double  quantum  well,  i.e.,  in  a  system  consisting  of  two 
different  quantum  wells  coupled  by  electron  tunneling.  We  will  concentrate  on  intersubband 
electronic  trzmsitions,  which  are  excited  by  far  ir  radiation  and  consider  both  optically  linear 
and  nonlinetir  effects.  We  will  focus  on  the  effect  of  the  light-induced  transfer  of  electrons  from 
one  well  to  the  other  one  in  the  double  quantum  well.  It  is  important  that  in  biased  quantum 
wells  such  a  transfer  can  occur  against  the  electric  field  force  and  with  high  quantum  yield  (up 
to  0.5). 


A  quantum  well  is  a  semiconductor  heterostructure  (see,  e.g.,  Ref.  1)  whose  potential 
confines  electrons  to  a  small  region.  Such  confinement  brings  about  quantum  splitting  of  the 
electron  energy  bands  into  subbands  separated  by  excitation  energy  on  the  order  of  /m* a?, 
where  m*  is  the  electron  effective  mass  and  a  is  a  confinement  size  (width  of  the  well).  In  what 
follows,  we  will  assmne  that  the  conduction  band  states  2ire  populated  due  to  a  modulation 
doping  of  the  barrier  regions  and/or  an  incoherent  optical  excitation  from  the  valence  band, 
and  consider  purely  electronic  transitions  between  subbands  of  the  conduction  band  [often 
called  QWEST  (Quantum  Well  Electronic  (inter )Subband  Treinsitions)].  We  will  also  assume 
the  electron  density  to  be  small  enough  to  exclude  excitonic  amd  other  mamy-body  effects. 


Much  work  hats  been  done  on  the  electronic,  optical  aind  kinetic  properties  of  semicon¬ 
ductor  double  quantum  wells  (see,  e.g.,  Ref.  1  and  references  cited  therein,  amd  adso  recently 
published  papers^“*^,  which  aie  relevant  for  the  present  work).  A  fundamental  phenomenon, 
which  is  a  subject  of  the  study,  is  resonamt  tunneling  between  the  quamtum  wells.  A  distinc¬ 
tive  feature  of  this  phenomenon  is  a  considerable  enhamcement  of  the  tunneling  probability 
if  the  energies  of  the  donor  and  acceptor  levels  are  close  enough.  To  describe  theoreticadly 
this  phenomenon  or  interpret  experimented  results,  most  of  the  above  cited  works  rely  on  the 
use  of  the  Schrodinger  equation.  In  this  approach^,  the  wave  functions  of  resonant  levels  in 
the  wells  aire  mixed  due  to  tunneling,  amd  these  states  repulse  forming  a  doublet  sepaurated 
by  the  energy  2|r|,  where  r  is  the  tunneling  aunplitude.  The  tunneling  is  described  by  the 
delocalization  of  the  electron  wave  function.  Such  tunneling  is  often  called  coherent,  amd  we 
will  follow  this  terminology.  It  is  well  understood  (see,  e.g.,  Refs.  3,  4  amd  7)  that  relaucation 
destroys  coherence,  and  madces  tunneling  incoherent  (stepwise).  When  the  relaxation  rate  T  odat 
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becomes  on  order  of  the  tunneling  amplitude  r  or  greater,  the  incoherent  (stepwise)  tunneling 
takes  place  between  non-mixed  states. 

To  describe  a  general  case  of  an  arbitrary  relaxation  we  will  use  the  density  matrix 
technique,  which  allows  one  to  fully  take  into  account  the  relzocation,  including  the  dephasing 
contribution  to  the  polarization  relaxation  rate.  Such  contribution,  which  is  usually  neglected, 
may  be  important.  The  Schrodinger  equation  approach  will  also  be  used  below  to  describe 
coherent  tunneling  for  the  sake  of  comparison. 

To  explain  the  essence  of  the  electron  transfer  elfect,  let  us  consider  an  asymmetric  double 
quantum  well  with  an  electric  field  applied  perpendicular  to  the  well  plane.  The  schematic  of 
the  confining  potential  and  electron  levels  (subbands)  is  shown  in  Fig.  la  with  |1)  and  |2)  as 
the  ground  states  in  the  narrow  (N)  and  wide  (W)  wells,  respectively.  The  excited  state  in  the 
N  well  is  |3),  and  in  the  W  well  |4).  Let  us  assume  that  the  electric  field  aligns  the  excited  levels 
|3)  and  |4),  so  that  tunneling  from  one  excited  level  to  the  other  one  is  resonantly  enhanced. 

Qualitatively,  the  electron  transfer  effect  is  most  pronounced  in  the  coherent  tunneling 
case.  We  should  mention  that  there  exists  a  convincing  evidence  of  feasibility  of  at  least 
partially  coherent  tunneling,  based  on  the  observation®  of  coherent  oscillations  in  a  double 
quantum  well,  which  has  been  suggested  earlier®.  Also,  the  possibility  of  partially  coherent 
timneling  is  witnessed  by  the  observations  of  the  electron  tramsfer  kinetics  independent  of  the 
barrier  width®,  and  of  the  resistance  resonaince  in  coupled  quamtum  wells'®.  However,  it  is 
not  known  how  strongly  the  observed  effects  aire  changed  if  relatxation,  especiaiUy  polarization 
relaxation,  is  appreciable.  Below  we  shadl  describe  a  coherent  pictmre  of  the  electron  tramsfer 
effect  and  address  a  general  caise  in  Sec.  3. 

In  the  coherent  caise,  the  aligned  excited  states  form  a  doublet,  the  upper  amd  lower 
components  of  which  we  denote  as  |4-)  amd  |— ).  The  |±)-state  wave  functions  are  delocalized 
over  both  the  N  amd  W  wells  due  to  resonant  tunneling.  In  contraist,  the  lower  levels  are  not 
adigned,  amd  the  |1)  state  is  baisically  localized  in  the  N  well  amd  |2)  in  the  W  well.  Since  the 
subbamd  splitting  of  the  W  well  is  smadler,  the  overall  ground  state  is  |1)  in  the  N  well  (see 
Fig.  lb).  We  aissume  both  the  electron  density  amd  temperature  to  be  not  very  high,  so  that 
only  the  |1)  state  is  considerably  populated. 

Suppose  that  ir  light  excites  an  intersubband  tramsition  in  the  N  band,  i.e.  one  of  the 
transitions  of  the  type  |1)  — ♦  1±)  shown  in  Fig.  lb  by  a  wavy  arrow.  Since  the  splitting  of 
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the  levels  in  the  N  well  is  eisstimed  to  be  considerably  greater,  the  radiation  does  not  excite 
a  transition  in  the  W  well.  The  electron  excited  to  either  of  the  |±)  states  is  quzmtum- 
mechanically  delocalized  over  both  the  wells.  Subsequent  relaxation  brings  about  electron 
transitions  to  the  ground  states  |1)  and  j2)  shown  in  Fig.  lb  by  dashed  arrows.  The  transition 
rates  are  proportional  to  the  probabilities  for  an  electron  to  be  localized  in  the  corresponding 
wells  and,  for  aligned  levels,  are  on  the  same  order  of  magnitude.  Thus,  with  an  appreciable 
probability,  the  electron  comes  to  the  state  |2),  which  is  mainly  locadized  in  the  W  well. 

Summarizing,  a  net  result  of  the  photoexcitation  of  the  intersubband  transition  in  the  N 
well  is  a  transfer  of  the  electron  from  the  N  well  to  the  W  well  in  the  direction  of  the  potential 
increase  (see  Fig.  lb),  i.e.  against  the  direction  of  of  the  field  force.  Indeed,  the  energy  needed 
for  such  a  transfer  is  talcen  from  the  exciting  radiation.  Note  that  if  the  transition  in  the  wide 
well  is  excited,  the  electron  transfer  would  occur  in  the  direction  of  the  field  force. 

The  closest  counterpart  of  the  above  described  effect  is  the  observation  by  Sauer,  Thonke 
and  Tsang^  of  photoinduced  space-charge  buildup  due  to  asymmetric  electron  and  hole  tunnel¬ 
ing  in  coupled  quantum  wells.  The  effect  of  Ref.  2  is  similar  to  the  present  effect  in  regard  to 
electron  transfer  against  the  electric-field  force  but,  nevertheless,  is  essentially  different  in  the 
following  respects.  First,  there  is  no  relaxation  involved  in  charge  buildup  in  Ref.  2,  and,  as  a 
result,  the  electron  buildup  is  minimum  for  the  levels  aligned,  while  in  our  case  it  is  maximum. 
Also,  for  the  aligned  excited  levels  after  switching  off  the  optical  excitation,  the  charge,  which 
has  been  transferred  between  wells,  disappears  in  a  time  on  the  order  of  the  resonant  tunneling 
time,  while  in  our  case  the  charge  transferred  is  stable  on  this  temporal  scsde.  Second,  the 
effect^  is  induced  by  interband  tr2msitions,  and,  therefore,  the  portion  of  the  photon  energy 
accumulated  in  the  potential  energy  of  a  transferred  electron  is  smadl,  as  distinct  from  the 
present  effect  based  on  intersubband  transitions.  Third,  the  charge  trzmsfer  in  Ref.  2  is  based 
upon  the  difference  in  the  tunneling  time  of  the  electrons  in  the  conduction  band  and  holes 
in  the  valence  band,  while  no  conduction-band  holes  participate  and  no  such  requirement  is 
relevant  for  the  present  effect. 

The  Schrodinger  equation  description  is  presented  in  Sec.  2.  The  density  matrix  approach 
is  developed  in  Sec.  3,  including  obtaining  the  basic  equations  (Sec.  3.1),  study  of  the  temporal 
dynamics  (Sec.  3.2),  stationary  solutions  (Sec.  3.3)  and  numerical  examples  (Sec.  3.4).  The 
results  obtained  in  the  paper  are  discussed  in  Sec.  4. 
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2.  COHERENT  ELECTRON  TRANSFER 
IN  SCHRODINGER  EQUATION  FORMALISM 

An  electron  in  a  heterostructure  is  characterized  by  the  quantum  number  of  the  state  in 
the  well  i  =  1,2,+,—  and  the  momentum  p  of  the  movement  in  the  well  plain,  its  energy  being 
Ci  +  p^/2m*,  with  Ci  as  the  subband  edge  energy.  For  typical  times  of  interwell  tunneling,  which 
are  normally  much  less  then  the  electron  translation  relaxation  times,  one  can  consider  p  as  a 
conserving  quantum  number.  Also,  the  photon  momentum  is  much  less  then  a  characteristic 
electron  momentum,  which  allows  one  to  consider  p  as  conserved  by  the  interaction  with  light. 
Teiking  into  account  that  subb2mds  in  the  well  tire  highly  parallel,  one  concludes  that  the 
energy  of  an  intersubband  transition  j,  p  — » i,  p,  either  tunneling  or  electromagnetic,  does  not 
depend  on  p,  and  is  simply  the  transition  energy  for  the  states  in  the  well,  eij  =  e,  —  Cj.  These 
arguments  allow  one  to  describe  electron  states  in  the  well  eind  kinetics  of  the  intersubband 
transitions  separately  from  the  lateral  movement,  as  conventionally  done. 

The  maximum  counter-field  transfer  effect  occurs  at  low  temperatures  T  <C 
which  we  assume.  The  rate  equations  describing  the  populations  n,  of  the  states  |i),  i  = 
1,2,—,+  have  the  form 


dni 

'W 


-to±ni+  (7i±  +  w±)  n±  +  7i2n2 


dn2 

=  72±«±  -  7l2»l2 


dn± 


=  u;±ni  -  (u>±  +  7±)n±  . 


(1) 


Here  ±  means  either  +  or  ■jji  are  phenomenological  rate  constants  for  decay  from  higher- 
to  lower-lying  states,  |j)  — »  |j),  i  =  +,—,2,  j  =  2,1,  and  7±  =  72±  +  7i±;  and  w- 
are  well  known  Einstein  coefficients  for  the  transitions  |1)  — >  |+)  emd  |1)  |— ),  which  can 
conventionally  be  expressed  in  terms  of  the  corresponding  dipole  matrix  elements  and  d-i, 
w±  =  |Ed±i|^  r±i  +  jF^ij  ,  with  E  as  the  amplitude  of  the  light  wave  and  ft 

as  the  light  frequency.  To  simplify  notations,  we  use  the  system  of  units  in  which  h  =  1. 

The  probability  of  electron  transfer  from  the  N  to  W  well  is  equal  to  the  population 
number  n2.  Assuming  that  the  radiation  excites  one  of  the  transitions  |1)  —*  |+)  or  |1)  -♦  [— ), 
one  finds  from  the  stationary  solution  of  Eq.  (1)  that 


n2  =  72±w± 


727±  +  (272  +  72±) 


-1 


(2) 


where  72  =  712  is  the  decay  constant  of  the  state  |2). 
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For  the  limiting  case  of  nonsaturating  light  intensities,  w±  ^  727±(272  +  72±)~\  as¬ 
suming  both  trzmsitions  |1)  — *•  j-l-)  and  |1)  j— )  to  occur,  and  taking  into  accoimt  that  the 
contributions  of  these  transitions  axe  additive,  one  obtains  from  Eq.  (2) 

«2  =  72"*  (t"+72+/7+  +  .  (3) 


The  corresponding  expression  for  the  transfer  quantum  yield  Q  has  the  form 

Q  =  ^u;+72+/7+  +  u;_72_/7_)  -f 


(4) 


As  one  can  see  from  Eqs.  (l)-(4),  the  electron  transfer  kinetics  is  determined  by  the  decay 
constants  'jji  and  the  matrix  elements  d±i  which,  in  turn,  depend  on  the  wave  fvmction  mixing 
between  the  individual  wells.  This  mixing  is  described  by  the  probabilities  ’  and  P^  for 
an  electron  in  the  mixed  |±)  state  to  be  in  the  corresponding  N  or  W  well. 

To  find  the  decay  constants  7ij,  we  invoke  a  quantum- mechanical  idea  that  the  relaxation 
causes  localization,  and  an  electron  localizes  in  the  well  in  which  it  has  experienced  the  relax¬ 
ation.  This  assumption  is  valid  if  the  nonresonant  tunneling  rate  is  small,  which  is  equivalent 
to  neglect  of  the  overlap  of  the  wave  functions  in  different  wells  (see  also  below).  The  decay 
rate  of  the  excited  electron  in  the  N  well  is  equal  to  the  decay  rate  73  of  the  state  |3)  and, 
similarly,  the  decay  rate  in  the  W  well  is  74  of  the  |4)  state.  To  determine  d±i,  we  take  into 
accoxmt  that,  with  neglect  of  the  overlap,  the  electromagnetic  radiation  couples  the  groimd 
state  |1),  which  is  mainly  localized  in  the  N  well,  only  to  the  |3)  component  of  the  mixed 
states.  From  these  arguments  we  find 

7i±  =  734''’  .  =  ,  Id±,|^  =  |d„|'pf>  .  (5) 


An  important  characteristic  of  the  electron  transfer  process  is  the  mean  quantum  yield 
Q  for  the  wide-band  radiation,  i.e.  for  light  whose  spectral  width  is  much  greater  than  the 
radiation  transition  widths.  It  is  defined  by  the  expression  [cf.  Eq.  (4)] 


Q 


=  72^  n2dfi  J 


^  -I 


{w+  +  w-)dQ, 


(6) 


From  this,  taking  into  account  Eqs.  (3)  and  (5),  one  obtains 


Q  =  74  Pf’Pr’/7+  +  J’l"’Pi"77- 


(N)  d(W) 


,(N) 


+  Pr 


(N) 


1-1 


(7) 
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To  determine  the  transfer  probability  (2),  we  need  also  to  estimate  the  interwell  transition 
constant  72  =  712.  Note  that  72  is  proportional  to  the  probability  Pf  ^  electron  in  the 

state  |2),  which  is  mainly  loccilized  in  the  W  well,  to  appear  in  the  N  well  and  experience  a 
relaxation  there.  Assuming  the  relaxation  of  all  states  in  the  N  well  to  occur  with  the  same 
rate,  determined,  e.g.  by  collisions,  one  can  estimate  72  ~  73^2^^  ,  and  obtedn  from  Eq.  (2) 
the  satiurated  transfer  probability 

=  (1  +273P2^''V74Pr^)~‘  .  (8) 

Since  the  nonresonant  transfer  rate  is  much  less  than  the  resonant  one.  i.e.  Pj^^  «pr\the 
value  of  is  close  to  unity. 

We  have  numerically  solved  the  stationary  Schrodinger  equation  in  the  coordinate  repre¬ 
sentation  with  the  confining  potential  V{x)  shown  in  Fig.  1.  For  the  method  of  solution,  the 
dependence  of  V’(i)  on  the  chemical  composition  of  the  well  eind  effective  mass  m*  adopted  in 
the  calculations,  see  Ref.  5. 

As  an  example,  we  consider  the  double  well  system  (see  Fig.  1)  AW  barrier/W  well/WN 
barrier/N  well/NB  barrier  with  the  widths  equal  to  100/19/8/14.5/100  nm,  all  the  barriers 
being  Alo.1Gao.9As,  which  corresponds  to  the  well  depth  Uq  =  77.5  meV.  The  solution  yields 
the  dependence  of  the  excited  level  detuning  €43  =  2.56(E  —  6.4)  meV,  where  E  is  the  bias  field 
in  kV/cm;  the  aligning  electric  field  is  Eq  ~  6.4  kV/cm.  For  the  aligning  field,  the  solution 
gives  the  energy  levels  (meV)  ci  =  4.7,  €2  =  15.9,  e_  =  40.9  and  e+  =  43.9;  the  tunneling 
amplitude  r  =  1.5  meV;  the  dipole  matrix  elements  d+i/e  =  2.4  nm,  d_i/e  =  —3.0  nm  with 
e  as  the  elementary  charge;  and  the  localization  probabilities  Pj^^  =  0.98,  P^^^  =  0.99, 
P^^^  =  0.58  and  pj.^^  =  0.42.  From  the  last  set  of  data  we  see  that  the  state  |1)  is,  in  fact, 
localized  in  the  N  well  and  |2)  in  the  W  well,  while  the  |-|-)  aind  ]— )  states  are  almost  evenly 
delocalized  over  both  the  wells,  as  assumed  above. 

The  transition  rate  72  with  P2^^  =  0.01  (see  above)  is  very  small  with  respect  to  73, 
which  ensures  a  high  saturated  probability  (8)  w  0.99,  low  optical  excitation  rates 
needed  to  achieve  saturation  of  the  |2)  state,  ~  72  •C  7±,  and  comparatively  long  Hfetime 
tc  =  7^^  of  the  transferred  charge  after  switching  off  the  radiation.  In  practical  terms,  the 
typiczd  decay  rate  of  the  excited  states  is  7  =  10^^  s“^  =  0.66  meV,  which  yields  tc  =  0.1  ns. 
For  the  linewidth  of  ~  2  meV  characteristic  of  QWEST  and  the  dipole  elements  given  above, 
the  saturation  light  intensity  =  e±i72/o’±  can  be  estimated  as  ~  60  kW/cm^. 
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Besides  the  data  shown  above,  the  computation  provides  escape  rates  "ysi  from  the  |i) 
states  in  the  double  well  to  the  B  region  (Fig.  lb),  i.e.  in  the  direction  of  the  potential  drop.  As 
expected,  the  largest  of  the  obtained  rates  are  those  for  the  excited  states,  7b_  =  1.8  ^eV  and 
7B+  =  1.1  /xeV.  These  escape  rates  play  the  role  of  the  rate  constants  for  a  parasitic  process 
of  the  light-induced  leak  from  the  quantum  well.  However,  comparing  7b±  to  the  tunneling 
amplitude  r  =  1.5  meV  and  also  to  the  rates  72±  ~  0.3  meV  of  the  population  of  the  (2)  state 
from  |±),  we  arrive  at  the  conclusion  that  the  escape  current  is  negligibly  small  with  respect  to 
the  interwell  tunneling  current,  and  can  cause  only  a  very  small  positive  charging  of  the  well 
system  as  a  whole  without  affecting  the  counter-field  electron  transfer. 

The  mean  transfer  quantum  yield  Q  as  a  function  of  the  bias  electric  field  E  calctilated 
from  Eqs.  (5)  and  (7)  for  73  =  J4  is  shown  in  Fig.  2  by  a  solid  line.  As  one  can  see,  Q  has  a 
rather  sharp  resonance  at  the  field  E  =  6.4  kV/cm,  which  exactly  corresponds  to  alignment 
of  the  excited  levels  in  the  two  coupled  wells.  The  maximum  value  is  Qmax  ^  0.55  and,  as 
the  computations  show,  it  essentially  does  not  depend  upon  the  barrier  width,  which  is  a 
consequence  of  the  coherency  of  the  tunneling.  As  one  cam  conclude  from  Eq.  (7),  the  most 
favorable  case  for  the  counter-field  transfer  is  not  the  simplest  choice  73  =  74  as  above,  but 
rather  74  >73,  in  which  caise  Q  cam  be  close  to  unity. 

The  tunneling  amplitude  r  is  an  essential  pairameter,  which  determines  whether  the  tim- 
neling  is  coherent  or  not  (see  Secs.  3  and  4).  This  amplitude,  which  is  equal  to  half  the  minimum 
value  of  the  doublet  splitting  c+_,  has  been  determined  from  the  solution  of  the  Schrodinger 
equation  ais  a  function  of  the  width  Lwn  of  the  interwell  (WN)  barrier  for  the  system  of  the 
above  type  with  the  widths  100/19/£w';v/14.5/100  nm.  The  result  is  shown  in  Fig.  3,  from 
which  we  cam  see  that,  in  the  range  of  Lwn  considered,  r  is  am  exponential  function  of  LwN^ 
as  expected. 

The  numericad  solution  given  above  illustrates  baisic  features  of  the  counter-field  transfer 
effects,  amd,  in  principle,  takes  into  account  mixing  between  all  the  states  in  both  the  wells, 
amd  effect  of  the  continuum  states.  However,  the  nature  of  the  numerical  solution  is  such 
that  the  amailyticad  dependence  of  the  effect  on  the  pauraimeters  of  the  problem  remain  uncleair. 
To  elucidate  this  dependence,  we  will  employ  a  conventional  (see,  e.g..  Ref.  1)  approximate 
approaich  baised  on  the  mixing  of  only  two  excited  states,  |3)  amd  |4)  amd  neglect  of  the  overlap 
integral.  The  latter  is  equivalent  to  the  orthogonality  condition  (3|4)  =0.  The  transfer 
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integral  t  in  the  notation  of  Ref.  1  is  exactly  our  tunneling  amplitude  t,  and  the  approximation 
under  consideration  is  nothing  else  but  the  tight-binding  model  with  zero  overlap  in  the  basis 
of  the  two  states  |3)  and  j4}. 

The  stationary  Schrodinger  equation  in  this  basis  is  reduced  to 

(€3-e±)(3|±) +r‘{4|±)  =0  ,  (*4  -  «±)  {4|  ±  >  +  t  (3|  ±  )  =  0  ,  (9) 

with  as  eigenenergies.  The  solution  of  Eq.  (9)  has  the  familieir  form  e±  = 
^  64  -1-  €3  ±  ^€43  ,  where  £43  =  64  —  63  is  the  energy  mismatch  of  the  excited 

levels.  Also,  the  localization  probabilities  follow  from  Eq.  (9)  as 

pW  =  p(W)  =  1  _  p|W)  ^  J  _  pm  ^  ^ ^  ^  ^ -  (10) 

«43  +  +  4|Tp^  +4|r)^ 

Substituting  Eq.  (10)  into  (2)  we  find  the  analyticcd  expression  for  the  transfer  probability 

n2  =  \Gt\^  7472"'  (c^3  -f-  4t^)  (ft  -  6±i)^  +  I  (74C±3  -  73eT3)^  (11) 

Here  the  Rabi  frequency  G  =  Edia  and,  as  usual,  the  notation  tij  stands  for  Cj  —  Cj. 

The  mean  quantum  yield  of  transfer  is  found  from  Eqs.  (7)  and  (10) 

Q  =  (73  +74)74  [kl^  (73  +  74)^  +  7374C43  •  (12) 

As  we  see  from  Eq.  (12),  Q  as  a  function  of  the  level  mismatch  643  has  the  form  of  a  symmetric 
peak  with  width  kl  (74  +  73)  (7374)  maiximum  magnitude  Qmaz  =  (1  +  73/74)  •  For 

the  simplest  choice  73  =  74,  using  the  dependence  643  =  2.56(E  —  6.4)  meV  (see  above),  the 
mezm  quantum  yield  (12)  is  plotted  in  Fig.  2  as  a  function  of  E  by  the  dashed  line.  We  see 
that  both  in  magnitude  and  form  a  simple  analytical  formula  (12)  agrees  well  with  the  result 
of  the  complete  numerical  calculation  (Fig  2,  sohd  line).  This  shows  good  applicability  of  the 
truncated-basis  tight-binding  approach  used  above  to  derive  Eq.  (12). 

As  mentioned  above,  Qmax  =  0.5  for  73  =  74,  and  Qmaz  -*  1  for  74  >73.  The  values 
of  Qmaz  do  not  depend  on  the  tunneling  amplitude  r  nor,  consequently,  on  the  barrier  width. 
This  is  because  of  the  coherent  quanttim  tunneling  zissumed  above.  A  high  quantum  yield  of 
transfer  demonstrated  above  znd  its  dependence  on  the  relaxation  constant  ratio  deserves  a 
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physical  interpretation.  The  electron  transfer  counter  to  the  electric  field  force  is  based  on 
the  quantum-mechanical  delocalization  of  an  electron  over  both  the  wells.  This  delocalizatior 
counterplays  the  electric  field  force  and  dominates  in  the  caise  of  realistic  fields,  where  the 
above  approximation  is  valid.  As  a  result  of  the  quantum  delocalization,  the  electron  can  be 
foimd  in  the  W  well  with  a  50  percent  probability.  The  relaxation  in  the  W  well  with  rate  74 
causes  the  electron  to  get  localized  in  this  well,  i.e.  transferred  against  the  field  force.  The 
relaxation  in  the  N  well  with  the  rate  73  brings  the  electron  back  to  the  vround  state,  i.e.  is  a 
parasitic  process.  Therefore,  Q  «  1  for  74  73. 

3.  EXCITATION  AND  ELECTRON  TRANSFER 
IN  DENSITY  MATRIX  FORMALISM 

3.1.  Basic  equations 

In  the  previous  section  we  have  described  the  electron  excitation  kinetics  and  counter¬ 
field  transfer  using  the  Schrddinger  equation  formalism.  Such  an  approach  is  adequate  for 
the  coherent  mechanism  of  the  electron  tunneling.  The  main  drawback  of  this  approach  is 
the  neglect  of  the  polarization  relaxation  processes,  which  tend  to  destroy  the  coherence  of 
tunneling.  Below  we  present  a  general  theory  based  on  the  density  matrix  approach,  which 
allows  one  to  describe  the  full  range  of  the  interwell  tunneling  regimes  from  completely  coherent 
in  the  C2ise  of  small  polau'ization  relaxatioi*.  rate  to  the  opposite  case  of  completely  stepwise 
for  strong  polarization  relaxation.  The  simplifying  feature  of  our  approach  is  the  use  of  the 
tight-binding  model  in  the  restricted  basis  of  the  states  in  isolated  wells  |1) ,  |2) ,  |3) ,  (4)  (see 
Sec.  2)  and  of  the  releocation  constamt  model  for  the  relaxation  term  in  the  equation  of  motion 
for  the  density  matrix  (see  below). 

We  start  with  the  Hamiltonian  of  the  system  in  the  form  H  =  CiCtJcii  -I-  Vijolaj, 
where  aind  a  are  the  electron  creation  and  aimihilation  operators,  with  t,j  =  1,2, 3,4. 
The  one-electron  operator  V  describes  interaction  with  the  electromagnetic  field  and  electron 
interwell  tunneling,  and  its  independent  nonzero  elements  are  V31  =  — d3i  (Ec-'"‘-l-c.c.), 
V43  =  r.  The  one-electron  density  matrix  r  is  defined  ais  Its  diagonal  matrix 

elements  are  the  population  probabilities  n,-  =  r,-,-. 

The  equation  of  motion  for  r  can  be  obtained  in  the  ususd  way  by  commuting  the  pair 
operator  aja,-  with  the  Hamiltonian  and  adding  the  relaxation  term.  This  has  the  well-known 
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form 


|^=t[r,c  +  F]-^  .  (13) 

where  the  one-electron  energy  operator  e  is  defined  els  {i  je]  j)  =  and  R  is  the  relaxation 

operator.  In  the  low-temperature  case,  i.e.  for  neglect  of  thermal  activation,  the  diagonal  part 
of  R  describes  spontcineous  decays  from  higher-  to  lower-lying  levels,  and  in  the  the  model  of 
relaxation  constants  has  the  form 

Rfi  —  ^  ^  ^  ^  '  (1^) 

i  (j<0  i  (>>•) 

where  is  the  rate  constEint  for  spontaneous  decay  |ji)  — +  |i). 

In  what  follows,  we  will  neglect  direct  relaxation  transitions  which  involve  nonresonEint 
interwell  tunneling,  |3)  |2)  Eind  |4)  — >  |1),  on  the  ground  of  small  probability  of  nonresonsmt 

tunneling  with  respect  to  the  resonant  one.  The  rates  of  the  above  processes  axe  negligible 
with  respect  to  the  rates  of  the  collateral  two-step  processes  involving  the  resonant  tunneling, 
13}  |4)  — ►  |2)  and  |4)  — >  13)  — ♦  ll),  which  will  be  taken  into  account.  Such  a  neglect  is 

equivalent  to  the  Eissumption  of  IocelI  character  of  the  relaxation  [see  the  discussion  preceding 
Eq.  (5)|.  However,  the  nonresonant  tunneling  process  12)  ll)  should  be  included  despite  its 
small  rate,  because  there  is  no  resonant  process  to  compete  with  it.  Thus,  the  following  decay 
rate  constants  should  only  be  taken  into  consideration:  73  =  713,  74  =  724,  72  =  7i2*  These 
rate  constants  have  the  sEime  meaning  els  in  Eq.  (5). 

The  nondiagonal  part  of  R  describes  the  polarization  relaxation  and  in  the  model  under 
consideration  is  given  by 

R-ij  —  1  ~  2  ’  (13) 

where  >  0  is  the  pure  dephasing  term. 

Let  us  introduce,  in  the  usual  way,  the  density  matrix  in  the  interaction  representation 
p  =  exp(— zet)r  exp(tet).  The  equation  of  motion  for  p  follows  from  Eq.  (13), 

^=.1p,c?1-s  .  (16) 

Here  U  =  exp(tct)V  exp(— tct)  and  5  =  exp(tct)ilexp(— tc<)  Eire  the  interaction  Emd  reltixation 
operators  in  the  interaction  representation.  The  diagontd  part  of  5  obviously  coincides  with 
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that  of  R  as  given  by  Eq.  (14).  The  nondiagonad  pcirt  of  R  is  simileir  to  Eq.  (15)  with  the 
substitution  pij  for 

The  equations  for  the  density  matrix  elements  follow  from  Eq.  (16)  taking  Eqs.  (14)-(15) 
into  account: 


dni 

dt 

dn2 

W 

dnz 

dt 

dn4 

'dT 


=  -2Im[/t)i3V3i  exp(iC3it)]  +  72712  +  73^3  , 

=  -72W2  +  74^4  , 

=  2Im[px3V'3i  exp(te3it)]  -73713  +  2Im  [p43r*  exp  (-7^43*)]  , 
= -2Im[p43r*  exp  (-7643^)] -74714  , 


^Pl3 

dt 

dpu 

dt 

dp43 

dt 


=  7(711  -  ^73)  VVi  exp  (-7631*)  +  7pi4r  exp  (7643*)  -  ri3Pi3 
=  -*P43^3*iexp(“*^3lO  +  */>13r*  exp  (-7643*)  -ri4/9l4  , 
=  *Pl4^3*  exp  (-7631*)  +  t  (n4  -  713)  T  eXp  (76430  “  r43*>43 


(17) 


The  system  of  equations  (17)  is  exact  for  the  model  under  consideration.  Below  we  adopt 
the  resonant  approximation,  which  is  also  synonymously  called  the  Rotating- Wave  Approx¬ 
imation  (RWA).  Applicability  of  this  approximation  is  well  established  for  optical  fields  not 
very  strong.  Techmcally,  RWA  is  equivalent  to  taking  into  accoimt  only  the  terms  containing 
slowly-oscillating  temporal  exponentials  in  Eq.  (17).  Doing  so,  we  can  explicitly  determine  the 
temporal  dependence  of  the  polarizations, 


Pl3  =  Pl3exp[7(fi  -  631)*]  ,  Pl4  =  ^»14exp(7(n  -  641)*]  ,  P43  =  P43  exp  [7643*] 


(18) 


and  find  the  system  of  equations  with  constant  coefficients 


dni  /  - 

—  — 2Im(pi3G)  -I-  72772  +  73^3  , 
dti 

=  2Im(  P13G)  -  73773  +  2Im(p43r*) 

(tIi  —  773) -f- 7rpi4  —  gi3Pi3  , 

dp43 

dt 


—  IT  (774  —  773)  -f  lG*p*^  —  5743^43  , 


dn2 

~dt 

dn4 

~dt 

dpi 


=  -72”2  +  74*34  , 

= -2Im(p43r*)  -  74774  , 


dt 


-  =  ~iG*pl3  +  iT*p\z  —  g\4P\\ 


where  the  notations  are  introduced  [G  is  defined  after  Eq.  (11)] 

fifl3  =  Tis  -|-7(n  -  631)  ,  =  Tu  +  *(ft  -  641)  ,  g43=r43+te43  • 


(19) 


(20) 
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3.2.  Temporal  dynamics 

The  fundamental  system  of  equations  (19)  describes  both  the  temporal  evolution  amd  sta¬ 
tionary  levels  of  the  electron  polairization  and  population  numbers.  Beginning  with  dynamics, 
the  simplest  effect  is  known  to  be  the  quantum  beats:  the  system  is  excited  with  a  short  pulse 
of  light  amd,  after  the  pulse  is  over,  the  population  numbers  and  polarizations  are  changing 
in  time,  in  some  cases  in  an  oscillating  manner.  Such  observation  conditions  are  in  a  close 
correspondence  with  the  experimentad  study  of  Ref.  8. 

The  kinetics  of  the  quantum  beats  is  described  by  Eq.  (19),  where  G  =  0.  In  this  case, 
we  see  that  the  equations  related  to  the  tunneling-connected  states  |3)  and  |4}  decouple  from 
the  rest  of  equations  (19),  forming  a  closed  system 


dn^ 

dt 


-73^3  +  2Im(r*/>43) 


drii 

dt 


— Un4  -  2lm(r*/J43) 


dp43 

dt 


=  IT  (n4  —  Us)  —  943P43  • 


(21) 


The  solution  of  this  system  is  given  by  a  superposition  of  exponentials  with  the  coefficients, 
which  can  be  found  from  the  initial  conditions  determined  by  the  excitation  process.  However, 
the  corresponding  characteristic  exponents  An  are  the  eigenvalues  of  Eq.  (21),  which  do  not 
depend  on  the  initiad  conditions. 

The  system  (21)  is  equivalent  to  four  real  equations  (the  matrix  element  ^43  is  complex). 
The  corresponding  characteristic  equation  is  fourth-order  and  yields  analyticad  expressions 
for  the  four  eigenvadues  A„.  However,  these  expressions  are  too  complicated.  Therefore,  we 
will  assume  simplifying  conditions,  which  do  not  affect  the  generad  properties  of  the  solution. 
These  conditions  are  zero  level  mismatch  (C43  =  0)  and  equal  decay  rates  of  the  excited  levels 
(73  =  74)-  In  this  caise,  we  obtain  the  four  characteristic  exponents  Ai,  A2,  A+,  A_  to  be 


Ai  =  -73  , 


A?  =  — r43  > 


A± 


~2  +  ^43)  ± 


(22) 


We  cam  see  from  Eq.  (22)  that  the  first  two  exponents  correspond  to  purely  decaying  (ape- 
riodicad)  kinetics.  The  two  other  exponents  A±  may  have  am  imaginary  pairt  and,  consequently, 
describe  am  oscillating  evolution  depending  on  the  relation  between  the  pure  dephasing  term 
f43  [cf.  Eq,  (15)]  and  the  tunneling  amplitude  r.  Namely,  if  f43  <  4(r|,  then  the  oscillations 
of  the  populations  ns,  n4  amd  polarization  P43  with  the  frequency  wj  =  f4|r|^  —  jf43^ 
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are  possible,  depending  on  the  initial  conditions.  Let  us  emphasize  that  u;^  does  not  depend 
whatsoever  on  the  population  decay  rates  7,  . 

Physically,  the  oscillations  considered  above  mean  coherent  tunneling  transitions  of  an 
electron  from  one  quantum  well  to  the  other  and  back.  These  oscillations  decay  with  the 
rate  73  +  depending  on  both  the  population  decay  constant  73  and  the  dephasing  rate 

r43.  Such  coherent  oscillations  have  been  suggested  in  Ref.  7  and  observed  experimentally  in 
Ref.  8.  However,  in  the  theoretical  description  of  Refs.  7  and  8,  the  deph2ising  has  not  been 
taking  into  account.  If  the  dephasing  is  negligible  (f43  <C  |r(),  the  frequency  o/j  of  beats  is 
equed  to  the  doublet  splitting,  and  for  the  case  of  exactly  aligned  levels  we  have  Wfc  =  2  |t|,  in 
accord  with  Refs.  7  and  8.  However,  in  a  practical  situation,  the  dephasing  is  not  small,  and 
ujj,  <  2|r|.  As  the  results  obtained  above  show,  the  dephasing  relaxation  rate  r43  is  a  very 
important  parameter,  which  not  only  determines  the  values  of  the  frequency  and  decay  rate  of 
oscillations,  but  their  very  existence. 


3.3.  Stationary  solutions 

The  stationary  solution  of  Eq.  (19)  can  be  obtained  in  a  straightforward  manner.  The 
expressions  for  the  population  probabilities  of  the  excited  levels  are 


U4 


=  c(6  +  74)  -  0(0  +  74)  |c[6(3  + 74/72) +  274]  + 

7374  +  M73  +  74)  +  a  [73  (1  +  74/72)  -  a(3  +  74/72)  -  374]|  , 

=  o(73-o)  +  c6  |c [6(3 +  74/72)  +  274]  + 

7374  +  6(73  +  74) +  0(73(1 +  74/72) -a (3  + 74/72) -374]|  , 


(23) 


where  the  notations  are  introduced 

<.  =  2|Tp|GpRe(/-')  ,  6  =  2|7|^Re[(s,3Si4  +  |rnr'l  , 

c  H  2|GpRe  [(5^3,^14  +  |GH /"*]  ,  /  =  ff43  +  |7p) +5l3|Gp  . 

As  follows  from  from  the  second  of  the  equations  (19),  the  probability  n2  of  the  electron  transfer 
is  simply  related  to  n4  by 

n2  =  7472'^«4  •  (25) 
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The  expressions  (23)  are  exact.  However,  they  are  somewhat  cumbersome  and  difficult 
to  analyze.  In  particular,  it  is  difficult  to  see  from  Eq.  (24)  how  the  Schrodinger  approach 
results  (see  Sec.  2)  can  be  reproduced.  Therefore,  prior  to  giving  numerical  exaunples,  we  shall 
consider  simplifying  limiting  cases,  restricting  ourselves  to  an  opticcdly  nonsaturated  regime, 
kl-  Assuming  this,  we  obtain  from  Eq.  (23)  closed  expressions  for  the  population 

numbers 

”3  =  r'lsEu  —  (Q  —  e3l)(J^  —  f4l)]  [|T|^r43  (2ri4  —  74)  +74ri4  (r43  +  C43)]  + 

[Fis  (fi  —  641)  +  Fu  (Q  —  £31 )]  [|rp  (2F43  {Q  —  €41)  —  74C43)  +  74  (fl  —  f4i)  (F43  +  ^43)]  I  , 


2lrPlGI^  f 

”4  = - Y - [1^1^  +  Fi3ri4  —  (Q  —  €31)  (ft  —  ^41 )]  F43  (2Fi4  +  73)  + 

[Fi3  (ft  —  €41)  +  Fi4  (ft  —  €31)]  [73643  +  2F43  (ft  —  641)]  I  » 

(26) 

where  the  notation  is  used 

=  {  [kl^  +  Fi3Fi4  —  (ft  —  €31)  (ft  —  €41)]  +  [Fi3  (ft  —  641)  +  Fi4  (ft  —  C3l)r|  X 

[2|r|^F43  (73  +  74)  +  7374  (F43  +  643)]  . 

The  populations  (26)  are  proportional  to  \G^ ,  i.e.  to  the  light  intensity,  as  expected. 

From  Eq.  (26)  we  can  see  how  to  reproduce  the  results  of  Sec.  2.  If  the  polarization 
relaxation  rates  eire  much  smadler  than  the  tunneling  amplitude, 

Fi3»  Fi4,  F43  •<  |r|  ,  (28) 


the  denominator  X  (27)  in  Eq.  (26)  becomes  nearly  zero  for  ft  =  e+i  or  ft  =  e_i,  where  e±  are 
the  energies  of  the  doublet  levels.  In  this  case,  each  of  the  populations  n3  and  n4  as  a  function 
of  frequency  has  two  nonoverlapping  peaks  centered  at  the  frequencies  e±i .  In  the  vicinities  of 
these  peaks,  the  expression  for  734  hzis  the  form 


734  = 


2|rp|Gf  (Fi4e±3  —  Fi3C^3)  (73643  —  2F436:j:3) 


(4(rP  +  €^3)  (ft  —  e±i)^  +  (Fi46±3  —  Fi36:f3)^  [2|rpr43  (73  +  74)  +  C437374] 


(29) 


Taking  Eq.  (25)  into  account,  we  conclude  that  Eq.  (29)  agrees  with  the  resxilt  (11)  of  the 
Schrodinger  equation  approach  only  if,  in  addition  to  the  condition  (28),  pure  dephasing  is 
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negligible  with  respect  to  the  population  relaxation. 


ri4  74  1  Tu  «C  73  ,  r43  74  +  73  . 


(30) 


Inapplicability  of  the  Schrodinger  equation  approach  in  the  presence  of  the  dephasing  is  evident 
if  consider  the  mean  quantum  yield  Q  (6),  which  can  be  calculated  from  Eqs.  (25)  and  (29), 


Q  =  2  |r|^  r4374 


2r43  |r|^  (73  +  74)  +  7374^43 


(31) 


in  comparison  with  the  corresponding  result  (12).  Note  that  the  condition  (28)  may  be  realistic, 
but  (30)  is  normally  not  the  case  (see  below).  Therefore,  taking  into  account  the  poleirization 
relaxation  with  the  use  of  the  density  matrix  techniques,  is  essential. 


3.4.  Numerical  illustrations 

In  this  subsection  we  numerically  illustrate  properties  of  the  optical  excitation  in  the  cou¬ 
pled  quantum  wells,  with  emphasis  on  the  counter-field  electron  transfer.  We  aim  to  elucidate 
the  effect  of  the  polarization  relaxation  on  the  electron  excitation  and  transfer. 

The  exact  expressions  for  the  population  probabilities  733,  of  the  excited  levels  are 
given  by  Eq,  (23),  The  probability  n2  of  the  electron  transfer  in  accordance  with  Eq.  (25) 
simply  copies  734.  Since  the  expressions  (23)  crucially  depend  on  the  polarization  relaxation 
constants  ri3,  Fu,  r43,  we  begin  with  estimating  these  constzmts. 

Experimentally,  the  polarization  relaxation  constant  is  found  eis  the  width  of  the  optical 
absorption  line  under  nonsaturated  excitation  conditions.  Using  such  m  approach,  it  has  been 
established  in  Ref.  14  that  at  the  temperatures  T  <  100  K  and  the  transition  energy  e  <  100 
meV,  the  transition  width  F  does  not  significzmtly  depend  on  T,  and  F/e  w  0.03.  These  data 
show  that  the  contribution  of  the  optical  phonons  to  the  polarization  relaxation  under  the 
conditions  used  is  small.  The  optical  phonon  contribution  will  be  even  smaller  in  the  absence 
of  optical  phonon  emission,  which  is  the  case  for  C31  <  Wo,  where  Wo  is  the  optical  phonon 
frequency.  This  condition  is  satisfied  for  the  system  under  consideration. 

With  the  optical  phonons  excluded,  the  remaining  mechanism  of  the  dephasing  is  based 
on  the  fluctuation  of  plus-minus  one  monolayer  in  the  width  of  the  well.  This  mechanism  is  also 
supported^®  by  the  experimental  data  on  the  interband  transitions  in  quantum  wells.  Assuming 
this  mechanism,  we  arrive  at  a  simple  estimate  F/e  w  26a/ a,  where  a  is  the  well  width  and  6a  is 
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its  fluctuation,  6a  «  0.3  nm.  In  our  case,  e.g.  for  the  N  well,  a  =  14.4  nm  and  f  /e  ~  0.04.  For 
a  crude  estimate,  the  agreement  with  the  experimental  value  is  reasonable.  For  the  example 
considered,  this  estimate  yields  f  13,  fu  «  1.5  meV.  In  the  absence  of  optical  phonon  emission, 
a  typical  lifetime  of  the  intersubband  trcinsition  is  ~  3  ps.  This  yields  the  decay  constants 
73?  74  ~  0.2  meV.  Thus,  the  dephasing  contributions  to  the  polarization  relaxation  constants 
dominate,  in  contrast  to  the  applicability  condition  (30)  of  the  Schrodinger  equation  approach. 

To  estimate  the  interwell  polarization  relaxation  rate  r43,  we  notice  that  the  fluctuation 
of  the  excited  level  energies  £3  and  £4  are  always  greater  than  those  of  the  the  ground  states  £1 , 
£2.  In  this  case,  assuming  the  fluctuations  in  the  different  wells  to  be  independent,  we  obtain 
a  simple  relation  r43  =  (ri3  +  r|4)  ,  which  will  be  used  below. 

The  effect  of  the  polarization  relaxation  on  the  electron  excitation  can  be  traced  in  Fig.  4. 
We  can  see  in  Fig.  4a  that  in  the  case  of  a  weak  dephasing  (ri3  =  F^  =  0.1  meV’),  the  excitation 
contours  are  two  almost  separate  peaks,  which,  as  can  easily  be  verified,  are  positioned  at  the 
transitions  frequencies  t±  of  the  doublet  levels.  The  asymmetry  of  the  excitation  contour  is 
due  to  nonzero  level  mismatch  £43. 

We  emphasize  that  near  the  peak  maxima  114  ~  ns,  which  mezins  that  the  electron  transfer 
coimter  to  the  electric  field  force  occurs  with  a  high  probability,  em  electron  is  excited  by  the 
light  to  the  |3)  state,  but  appears  with  close  or  even  greater  probability  in  the  |4)  state  localized 
in  the  other  well.  Using  Eq.  (25),  we  can  make  sure  that  in  the  spectral  maxima  ns  «  1,  i.e.  the 
populations  are  saturated.  At  the  same  time,  the  parameter,  which  governs  the  polarization 
saturation,  |G|^  /F^j  <C  1.  This  means,  in  particular,  an  absence  of  field  broadening  and  low 
probability  of  excitation  into  the  continuum. 

With  an  increase  of  the  dephasing  rate  (see  Fig.  4b,c),  the  spectral  peaks  are  broadened 
and  overlap.  However,  the  population  number  n4  and,  consequently,  the  trsmsfer  probability 
n2  (25)  do  not  considerably  diminish.  This  witnesses  that  the  counter-field  transfer  effect 
persists  even  for  relatively  strong  dephasing.  For  Fi3  =  Fi4  =  10  meV,  the  doublet  structure 
is  completely  absent,  zind  the  absorption  contour  is  symmetrical  and  centered  at  the  frequency 
£31  of  the  transition  in  the  isolated  narrow  well  (Fig.  5c).  This  means  that  the  electron  is  first 
excited  within  the  N  well  and  then  tunnels  into  the  W  well,  i.e.  the  tunneling  is  incoherent. 

The  effect  of  the  polarization  relaxation  is  even  less  evident  for  the  case  of  nonequal 
dephasing  rates  in  the  two  wells  shown  in  Fig.  5.  With  an  increase  of  the  dephasing  in  the 
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wide  well,  the  peaks  begin  to  broaden  (Fig.  5a,  cf.  Fig.  4a),  but  then  this  broadening  stops 
(Fig.  5b)  and  reverses:  for  strong  relaxation  we  observe  a  single  ntirrow  peak  centered  at  the 
transition  frequency  631  of  an  isolated  narrow  well  (Fig.  5c).  Again,  this  fact  witnesses  that 
the  tunneling  is  incoherent.  The  narrowing  of  the  absorption  contour  for  the  case  of  strong 
dephasing  is  a  counterpart  of  the  w-.  ;]  known  spectroscopic  phenomenon  of  the  spectral  line 
collapse.  This  narrowing  can  be  understood  from  the  following  arguments.  With  zin  increase  of 
the  dephasing  in  the  W  well,  this  well  behaves  like  an  overdamped  resonator.  Such  a  resonator 
is  known  to  decouple  from  a  high-quality  resonator,  which  is  the  W  well.  In  such  a  way,  a  well 
with  a  strong  polarization  relaxation  does  not  considerably  perturb  the  other  well. 

To  focus  on  the  effect  of  counter-field  electron  transfer,  let  us  consider  Fig.  6,  in  which 
the  maximum  (in  the  light  frequency)  probability  of  this  transfer  is  shown  as  a  function  of 
the  tunneling  amplitude  r.  Note  that  for  the  double-well  system  under  consideration,  the 
magnitude  of  r  is  simply  related  to  the  thickness  Lwn  o{  the  interwell  barrier  (see  Fig.  3). 

In  the  case  of  low-to-intermediate  optical  saturation  (|G1  =  0.2),  we  see  from  Fig.  6a  that 
for  r  <C  Fia  the  transfer  probability  strongly  depend  on  t,  in  fact  n2  oc  |r|^.  This  a  feature 
of  a  noncoherent  electron  transfer.  With  increase  of  r,  the  probability  n2  levels  off.  The 
greater  713,  the  later  this  leveling  takes  place  and  the  lower  is  the  limiting  magnitude  of  n2. 
For  the  typical  value  r  w  1  meV,  the  transfer  probability  is  very  high,  n2  «  0.7.  Thus,  the 
counter-field  transfer  effect  is  strong  imder  the  realistic  conditions  considered.  For  a  high  light 
intensity  |G|  =  20  meV  (Fig.  6b),  the  transfer  is  completely  saturated,  n2  ^  1,  in  the  coherent 
tunneling  regime.  However,  there  exists  an  essential  difference  with  the  intermediate-intensity 
case  (Fig.  6a).  Namely,  the  effect  of  the  dephasing  is  opposite:  the  higher  Fia,  the  higher 
732,  and  the  sooner  the  saturation  sets  on.  This  counterintuitive  feature  stems  from  the  field 
broadening,  which  prevails  over  the  phase  relaxation  in  this  case. 

4.  DISCUSSION 

This  paper  pursues  two  interrelated  goals.  The  first  is  to  give  a  theoretical  description 
of  a  novel  effect  in  eisymmetric  double  quantum  wells,  namely,  the  electron  tramsfer  cotmter  to 
the  bias  field.  This  transfer  is  based  on  the  quamtum-mechanical  delocalization  of  the  electron 
over  the  resonant  states,  which,  in  the  caise  under  consideration,  prevauls  over  the  electric  field 
force  (note  that  the  light  itself  does  not  exert  any  significamt  force).  The  second  goal  of  the 
paper  is  to  develop  a  theory  which  describes  the  kinetics  of  the  electron  excitation  in  coupled 
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quantum  wells  taJcing  into  account  the  polarization  relzixation. 

An  approach  based  on  the  Schrodinger  equation  provides  a  comprehensive  description  of 
the  problem.  This  is  given  in  Sec.  2,  where  generaJ  expressions  (2)  for  the  electron  transfer 
probability  n2  and  (7)  for  the  mean  quantum  yield  Q  of  transfer  are  obtained  [see  also  Eq.  (12) 
and  Fig.  2]. 

However,  the  Schrodinger  equation  approach  is  only  valid  in  the  absence  of  the  polar¬ 
ization  relaxation  (dephasing),  which  always  exists  in  real  systems.  This  relaxation  destroys 
quantum-mechanical  coherence  and,  rigorously  speziking,  makes  the  stationary  states  of  the 
electron  nonexisting.  Physically,  the  polarization  releixation  originates  from  interaction  of  an 
electron  with  other  parts  of  the  system,  in  particular,  with  phonons  and  impurities  or  defects, 
which  are  not  taken  unto  account  in  the  Schrodinger  equation.  In  this  case,  the  density  matrix 
technique  is  adequate.  Such  technique  for  the  optical  excitation  of  a  double  quantum  well  is 
developed  in  Sec.  3. 

The  free  (i.e.  in  the  absence  of  the  optical  field)  time-dependent  solution  of  Eq.  (19)  [see 
also  Eq.  (20)]  gives  the  frequency  of  the  quantum  beats  ujj.  In  the  presence  of  phase  relaxation, 
this  frequency  does  not  coincide  with  the  spacing  of  the  excited  state  doublet,  and  is  determined 
by  the  dephasing  constant  f43  but  not  the  decay  constants.  If  the  phase  relaxation  is  strong 
enough,  r43  >  4|r|,  then  the  quantum  beats  disappear,  completely  replaced  by  a  monotonic 
relaxation.  The  results  discussed  may  be  of  interest  for  compeirison  with  recent  experiments  of 
Ref.  8,  in  which  quantum  beats  in  the  double  wells  have  been  observed. 

In  the  stationary  regime,  closed  analytical  expressions  (23)  and  (25)  for  the  populations 
n3  and  n4  of  the  excited  states  axe  obtained.  The  full  range  of  the  transition  from  a  nearly 
coherent  tuimeling  regime  to  an  almost  incoherent  one  can  be  traced  in  Figs.  4  and  5.  For 
a  wealc  polarization  relaxation,  the  excitation  profiles  reveal  a  two-peeik  structure  typicaJ  for 
the  coherent  tunneling.  In  this  case,  the  optical  wave  plays  the  role  of  a  probe  field  exciting 
the  system  to  the  doublet  states  |±).  As  the  dephasing  increases,  the  stationary  states  |±) 
are  no  longer  a  good  zero-order  approximation,  and  the  interwell  coupling  becomes  incoherent 
(stepwise):  the  first  step  is  the  excitation  from  the  state  |1)  to  |3)  of  the  W  well  followed  by 
the  second  step  of  the  tunneling  into  the  N  well.  In  qualitative  agreement  with  the  above 
picture,  with  an  increase  of  the  dephasing,  the  double-peak  structure  disappears,  replaced  by 
a  single  peak  centered  at  the  transition  frequency  C31.  If  the  dephasing  in  the  two  coupled 
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wells  is  increased  in  the  same  proportion,  the  excitation  profiles  are  broadened  (see  Fig.  4c). 
.\n  counter-intuitive  feature  appears  in  the  case  where  the  dephasing  is  increased  only  in  the 
acceptor  (W)  well:  eventually,  a  collapse  of  the  spectral  line  teikes  place  (see  Fig.  5c). 

To  avoid  possible  misunderstanding,  we  recall  that  the  form  of  the  excitation  profiles  in 
Figs.  4  and  5  is  obtained  under  the  assumption  that  only  the  transitions  from  the  groimd  state 
|1)  in  the  N  well  eire  excited.  These  transitions  are  positioned  on  the  energy  scale  to  the  right 
of  the  treinsitions  in  the  W  well.  Therefore,  in  practice,  the  right  parts  (fl  >  30  meV)  of  the 
excitation  profiles  should  be  close  to  the  ones  shown  in  Figs.  4  eind  5,  and  the  left  parts  may 
differ.  However,  for  the  nonsaturating  optical  excitation,  only  the  transitions  in  the  N  well  can 
occur,  because  only  the  |1)  state  in  this  well,  as  the  overall  ground  state,  is  populated  at  low 
temperatures.  Therefore,  under  these  conditions,  the  excitation  profiles  are  correctly  given  by 
Eq.  (26)  in  the  whole  spectral  region,  predicting  dependences  similar  to  those  of  Figs.  4  and  5. 

As  discussed  above,  the  phase  relaxation  drastically  affects  the  spectral  profiles  of  the 
optical  excitation.  Also,  it  is  understandable  a  priori  that  the  probability  n2  of  the  counter¬ 
field  electron  transfer  diminishes  with  the  increase  of  the  dephasing  because  the  transfer  is 
a  resonant  effect,  and  the  quality  of  the  resonance  is  reduced  by  the  polarization  relaxation. 
A  question  remains  to  which  degree  na  is  affected.  The  answer  is  contained  in  Fig.  6a:  for 
moderate  excitation  intensities  and  the  whole  conceivable  range  of  fi3  =  0,1  —  10  meV,  the 
probability  n2  changes  only  by  a  factor  of  0.3.  Therefore,  the  transfer  effect  itself  is  rather 
stable  with  respect  to  the  dephaising. 

Let  us  discuss  possible  experimented  observation  of  the  counter-field  electron  transfer. 
This  effect  can  be  detected  optically  by  monitoring  chemges  of  the  intersubband  absorption 
in  the  double  well:  as  the  transfer  proceeds,  the  absorption  band  shifts  to  a  lower  frequency 
by  the  amoimt  621-  Electrical  detection  of  the  treinsfer  is  also  possible.  In  this  case,  external 
conductors  should  be  in  contact  with  the  regions  A  and  B  in  Fig.  1.  However,  sichieving  the 
regime  of  a  stationary  current  in  the  external  circuit  is  problematic,  because  in  this  case  the 
barriers  AW  and  NB  should  be  penetrable  for  electrons.  If  so,  the  optical  excitation,  apart 
from  bringing  about  the  counter-field  electron  transfer,  would  also  increase  the  rate  of  the 
electron  escape  from  the  N  well  to  the  B  region,  i.e.  in  the  direction  favored  by  the  bias.  Thus 
the  coimter-field  transfer  may  be  completely  masked  by  this  leak  current. 

We  believe  that  the  most  reliable  is  the  detection  of  the  counter-field  transfer  based  on 
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the  capacitance  coupling  of  the  well  to  an  external  circuit.  Such  coupling  is  achievable  even 
with  the  thick  barriers  AW  and  NB,  thus  excluding  photoinduced  lezikage  from  the  N  well  to 
B  region  discussed  above. 

In  the  case  of  the  capacitance  coupling,  the  regions  A  and  B  (Fig.  1),  containing  a  dense 
electron  gas,  play  the  role  of  the  capacitor  plates.  For  the  regime  of  zero  current  in  the  external 
circuit,  the  counter-field  transfer  of  electrons  induced  by  switching  on  of  the  light  brings  about 
an  increase  of  the  potential  difference  AUab, 


^Uab  =  ^'Ktari'iAxI e  , 


Ax 


xdx 


(32) 


where  4',  is  the  wave  function  of  the  ji)  state,  e  is  the  mean  dielectric  consteint  of  the  well 
materizd,  a  is  two-dimensional  density  of  the  electron  gas  in  the  well,  and  x  hsis  meaning  of 
the  characteristics  distance  of  the  charge  transfer  (for  the  example  considered.  Ax  =  24  nm). 
We  emphzisize  that  AU  >  0  means  that  the  photocurrent  inside  the  well  is  directed  counter  to 
the  potential  drop. 

Alternatively,  if  the  capacitor  is  externally  kept  under  a  constant  potential  difference, 
then  an  exciting  light  pulse  brings  about  a  transient  current  in  the  external  circuit  opposite  to 
the  direction  favored  by  the  bias,  the  total  transferred  charge  being  Q  =  c5crn2Ax/L,  where 
L  is  the  AB  distance,  2ind  S  is  the  illuminated  area;  Ax/L  ss  0.1  for  the  system  considered. 

We  should  point  out  the  the  theory  presented  above  does  not  take  into  accoxmt  the 
photoinduced  electric  fields.  This  condition  can  always  be  met  if  the  density  cr  is  low  enough. 
In  this  case,  of  course,  the  potential  change  (32)  is  small.  However,  it  is  expedient  to  mention 
the  probable  quzditative  effects  of  the  photoinduced  fields.  The  potential  incrcEise  AUab-,  via 
changing  the  electric  field  inside  the  well,  affects  the  photoexcitation  and  electron  transfer. 
This  is  a  feedback  which  can  produce  enhanced  nonlineEu*  optical  responses,  similar  to  ones 
observed  in  Ref.  6  for  the  interband  transitions,  emd,  possibly,  zm  intrinsic  optical  bistability. 
We  shall  address  these  effects  elsewhere. 
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FIGURE  CAPTIONS 


Fig.  1.  Coupled  wide  (W)  and  narrow  (N)  quantum  wells  (a  -  in  the  absence  of  bias,  b 
-  in  the  presence,  where  the  excited  levels  are  aligned).  Schematic  of  the  confining  potential, 
energy  levels,  and  radiative  (wavy  arrow)  and  nonradiative  (dashed  arrows)  tremsitions.  The 
regions  A  emd  B  containing  a  dense  electron  gas  serve  as  electrodes  for  the  capacitance  couphng 
of  the  double  well  to  an  external  circuit.  The  insulating  barriers  AW  and  NB  are  supposed  to 
be  thick  and  high  enough  to  exclude  considerable  tunneling  through  them  (see  the  text). 

Fig.  2.  Qusmtum  yield  Q  of  the  electron  transfer  counter  to  the  field  force  as  a  function  of  the 
electric  field  E  applied  to  the  double  well.  The  data  are  obtained  from  numerical  computation 
according  to  Eq.  (7)  (solid  line)  and  with  the  use  of  analytical  formula  (12)  (dcished  line). 

Fig.  3.  Tunneling  amplitude  r  (meV)  as  a  function  of  the  central  barrier  width  Lwn  (nm) 
in  the  double  logarithmic  scale. 

Fig.  4.  Population  numbers  ns  (solid  line)  and  and  7x4  (dashed  line)  of  the  excited  levels 
as  functions  of  the  exciting  light  frequency  Q  (meV)  for  the  dephasing  constants  shown  in  the 
graphs.  The  dephasing  relaxations  in  the  two  coupled  wells  are  assumed  equed,  fis  =  ri4. 
The  other  relevant  parameters  are  |G|  =  0.2  meV,  73  =  74  =  0.2  meV,  r2  =  0.006  meV,  and 
€43  =  1  meV.  The  double- well  system  under  consideration  is  described  in  the  text. 

Fig.  5.  The  same  as  Fig.  4,  but  the  dephasing  in  the  narrow  well  is  the  same  for  all  the 
graphs  shown,  f  13  =const,  and  dephasing  in  the  wide  well  is  increasing,  f  u  =  0.1,  1,  10  meV. 

Fig.  6.  Treinsfer  probability  n2  calculated  in  its  spectral  maximum  as  a  function  of  the 
tunneling  simplitude  r  for  the  dephasing  rates  713  =  0.1,  1,  10  meV,  and  the  Rabi  frequencies 
G  =  0.2,  20  meV,  as  shown  on  the  graphs.  Other  parameters:  714  =  713,  the  decay  rates 
73  =  74  =  0.2  meV,  72  =  0.006  meV,  and  the  level  mismatch  £43  =  0. 
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